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Using 1/S expansion, we discuss the magnon spectrum of Heisenberg antiferromagnet (AF) on 
a simple cubic lattice with small dipolar interaction at small temperature T Tn , where Tn 
is the Neel temperature. Similar to 3D and 2D ferromagnets, quantum and thermal fluctuations 
renormalize greatly the bare gapless spectrum leading to a gap A ~ ojo, where uo is the characteristic 
dipolar energy. This gap is accompanied by anisotropic corrections to the free energy which make 
the cube edges easy directions for the staggered magnetization (dipolar anisotropy). In accordance 
with previous results, we find that dipolar forces split the magnon spectrum into two branches. This 
splitting makes possible two types of processes which lead to a considerable enhance of the damping 
compared to the Heisenberg AF: a magnon decay into two other magnons and a confluence of two 
magnons. It is found that magnons are well defined quasiparticles in quantum AF. We demonstrate 
however that a small fraction of long-wavelength magnons can be overdamped in AFs with S 1 
and in quantum AFs with a single-ion anisotropy competing with the dipolar anisotropy. Particular 
materials are pointed out which can be suitable for experimental observation of this long-wavelength 
magnons breakdown that contradicts expectation of the quasiparticle concept. 

PACS numbers: 75.10.Jm, 75.30.Ds 


I. INTRODUCTION 

The concept of elementary excitations (quasiparticles) is a powerful approach in the modern theory of many- 
body systemsAccording to this concept, each weakly excited state of a system can be represented as a set of 
weakly interacting quasiparticles carrying quanta of momentum k and energy fk- Processes of spontaneous decay of 
quasiparticles and interaction between them lead to a finiteness of their lifetime that is related to the quasiparticle 
damping Tk- It is reasonable to introduce the idea of quasiparticle only if its lifetime is sufficiently large or if the 
damping is much smaller than the energy (I\ <C Ck). As long-wavelength quasiparticles have the smallest energies, 
weakly excited states of a many-body system are represented as collections of long-wavelength elementary excitations. 
Thus, they should be well-defined according to the quasiparticles concept. As regards short-wavelength elementary 
excitations, they can be defined badly or even cannot exist at all for some momenta. This situation is realized, for 
instance, in liquid 4 He which has a termination point in its spectrum.—^ As short-wavelength elementary excitations 
are normally well-defined, particular systems with overdamped short-wavelength quasiparticles have attracted much 
attention in recent years4~— 

The quasiparticle concept is supported by many experiments in various systems and numerous microscopic calcula¬ 
tions in particular models. For example, it was found in Ref.d^ that Tk <C £k at k <C 1 and T <C Tn in 3D Heisenberg 
antiferromagnets (AFs) with a small single-ion anisotropy, where Tn is the Neel temperature. In particular, it was 
obtained that Tk ~ e|r 3 In r at S' ~ 1 and 1<Ct 3 , where r = T/Tn -C 1. For large spins S> 1, when the regime 
Tn/S <T< Tn exists, the damping is estimated as Tk ~ e^r 2 at k <C 1. 

On the other hand, it has been revealed recentl y 13 : 14 that small long-range dipolar interaction in 2D and 3D 
Heisenberg ferromagnets (FMs) makes a small fraction of long-wavelength magnons to be heavily damped^ It has 
been obtained that a peak appears in the ratio Tk/ek at very small but finite momentum. The peak height is of the 
order of unity even if the temperature is much smaller than the Curie one (i.e., when FM can be considered as a weakly 
excited system). This unexpected result contradicts the conventional wisdom about long-wavelength quasiparticles 
and expectation of the quasiparticle concept. It is demonstrated also in Refs d 3 : 14 that it is the long-range nature of the 
dipolar interaction that is responsible for the anomalous damping of some long-wavelength magnons. In the majority 
of real FM materials this effect is screened by magnetocrystalline anisotropy leading to the gap in the spectrum. 
Besides, it is difficult to observe the long-wavelength magnons breakdown experimentally due to very small values 
of the corresponding momenta. However recent progress in neutron spin-echo techniqu o l -d- holds out hope that the 
corresponding measurements will be carried out in suitable FM materials. 14 

The purpose of the present paper is to carry out similar analysis of the magnon spectrum in Heisenberg AF 
with dipolar interaction on a simple cubic lattice at T Tn using 1/S expansion. We obtain that similar to 2D 
and 3D FM s 13 : 17 quantum and thermal fluctuations lead to anisotropic corrections to the free energy which make 
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FIG. 1: Diagrams of the first order in 1/5 for self-energy parts discussed in the present paper. Diagram (a) comes from 
four-magnon terms © in the Hamiltonian whereas (b) stems from three-magnon terms ©. Lines in these diagrams stand for 
Green’s functions ©■ 


the cube edges easy directions for the staggered magnetization. These corrections to the free energy are naturally 
accompanied by appearance of a gap A in the magnon spectrum in the first order in 1/5. All these phenomena are 
of “order-by-disorder” origin. 

We obtain in accordance with previous results^— that dipolar forces split the magnon spectrum into two branches. 
Despite its smallness, this splitting is responsible for a considerable enhancement of the magnon damping: it opens a 
way for a decay of a magnon into two other spin waves and for a confluence of two magnons. These processes of decay 
and confluence contribute to magnon damping because the dipolar interaction leads to three-particle vertexes in the 
Hamiltonian. As a result, the main contribution to the damping arises in the first order in 1/5 from the loop diagram 
shown in Fig. ©b). This should be contrasted with Heisenberg AF, which does not have odd-particle vertexes due to 
the rotation invariance of the Heisenberg coupling. As a consequence, the magnon damping obtained in Ref^ 2 and 
mentioned above arises at finite T from four-magnon vertexes in the second order in 1/5. 

Our damping calculation in the first order in 1/5 shows that magnons are well defined quasiparticles if 5 ~ 1. 
In particular, we obtain a peak in the ratio Tk/ek at k ~ A /D, where D is the magnon velocity, which height is 
proportional to T/D <C 1. The peak height increases considerably at 5 > 1 in the regime Ty/5 -C T <C Tjv in 
which case Tk/fk ~ const and a fraction of long-wavelength magnons with k ~ A/D turns out to be overdamped. 
We demonstrate that the long-wavelength magnons breakdown arises also when a single-ion anisotropy appears in the 
system which competes with the anisotropy of the dipolar origin mentioned above. We argue that this phenomenon 
can be observed in cubic AFs TlMnF 3 and RbMnF 3 doped with a very small amount of cobalt. Notice that dipolar 
forces enhance greatly the magnon damping as compared with the purely Heisenberg AFs considered in RefJ 2 and 
mentioned above. 

The rest of the present paper is organized as follows. Sec. [TT] is devoted to Hamiltonian transformations and to 
description of the technique. Renormalization of the ground state energy and the real part of spectrum are discussed in 
Secs. mu and IIV1 respectively. The spin-wave damping is derived in Sec.fV] The long-wavelength magnons breakdown 
at 5 > 1 and in AFs with the competing single-ion anisotropy is considered in Sec. IVI1 Sec. IVIII contains our 
conclusion. Three appendixes are added with details of calculations. 


II. HAMILTONIAN TRANSFORMATION AND TECHNIQUE 
A. Hamiltonian transformation 


We discuss Heisenberg AF with dipolar interaction on a simple cubic lattice which Hamiltonian has the form 


l^m 


(i) 


where summation over repeated Greek letters is implied, Ji m = J > 0 for nearest neighbors and J; m = 0 for other 
couples of spins, 


n a/3 _ w 0 3 Rfm^lm ^apRf„ 

'°tIm 


4-7T 


f? 5 

n lrr 


( 2 ) 


is the dipolar tensor, 


Wq = 47T 


( g /^) 2 

Vo 


( 3 ) 


is the characteristic dipolar energy that is smaller than 1 K in the majority of magnetic materials, and vq is the unit 
cell volume. We assume in the present paper that wq <C J. By taking the Fourier transformation, one has from 
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Eq. © 


k 


( 4 ) 


where ,/ k = J2i Jim exp(tkR Zm ) and Q^ 3 = J2i Q?!l exp(ikR/ m ). 

It is convenient to represent spin components in the local coordinate frame as follows: S; = Sfx + (Sfy + 
Sfz) exp(fkoRi), where x, y , and z are mutually orthogonal unit vectors, exp(*koR;) describes AF spin ordering 
being equal to +1 and —1 on sites belonging to different magnetic sublattices, k 0 = (7r, 7r, n) is the AF vector, and 
we set the lattice spacing to be equal to unity. This representation allows one introducing only one sort of bosons via 
the Dyson-Maleev spin representation which has the form 


S? = 


sr = 


sr = s- 



a[ai 


(5) 


Using the Holstein-Primakoff transformation instead of the Dyson-Maleev one © does not change the results obtained 
below. Taking the Fourier transformation in Eqs. © and using the relation S k = S^x + S^. +ko y + S^ +ko z one obtains 

from Eq. © for the Hamiltonian H = E gs + XEi where 


E qs = ~z:S 2 J 0 N 


( 6 ) 


is the classical ground-state energy, N is the number of spins in the lattice, and Hi denote terms containing products 
of i operators a) and a. In particular, Hi = 0 because it contains only Q^ and with a ^ (3 which are equal to 
zero. One has for the rest terms which are essential for further consideration 


H 2 

Hz 

Hi 


E ^ k « k «k + ^ (a k «- k 

'H £ (•«:«. 

ki+k 2 +k 3 =0 


t t \ „ t f? k B* t t \ 

+ a k a -k ) + £ka k+ko a k + — a k a_ k+ko + “Tf a k a - k+ko I > 

a kl ( a L _ a ~k 2 ) a - k3 + Q k2 a kl+ko (a k2 + a_ k2 ) a_ k3 ) , 


4/V E ( a k i a -k2 a -k 3 a - k4 (~ 2J k 2 + Qk 2 - Qk 2 +k 0 ) 

ki+k 2 + k3 + k4 =0 

+ a kl a k2 a -k 3 a -k 4 (-2J kl+k3 + Qki + Qk i+ ko ~ 2 Q k I+ k3 + ko ) 

- < + k 0 ((QZ - QZ+J < + «- k2 [QZ + QZ + J) «- k3 «- k4 ), 


(7) 

( 8 ) 

(9) 


where 

Ek = sj 0 -|(Qr+of +ko ). 
Bk = sjk - 1 (or - Q k + ko ), 
fk = (QZ + k 0 ~Q x k v ), 

Bk = if (QZ + k 0 + QZ) ■ 


( 10 ) 
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B. Green’s functions 


It is convenient for further calculations to introduce the following set of retarded Green’s functions: 


G(w, k) = (a k , a k ) w , 

F(cu, k) = (a k ,a-k) w , 
Q(uj,k) (a k + ko , a k )tj > 
^*(cJ,k) — (^k+ko 5 k)w ? 


G(u,k) = <al k ,a_ k }„ = G*(-w,-k), 
^(w.k) = (aL k ,a k ) w = F*(-u, -k), 

G(u, k ) = (a^ k+ko ,a_ k } w = k), 

J-t( w ,k) = (fll k+kol 4) u =^(-w-k). 


We have two sets of Dyson equations for them one of which has the following form: 


( 11 ) 


^E k — to A -Ek Bk A n k — 8k + 5 k — Bk A F k ^ 

Bk A II k E k + ui A Ek Bk A Vl 8k + 5 k 


( g\ 

Ft 


fo 1 ) 

<t k T 5 k +ko Bk T T’k-fko E k +ko ^ A Tk-|-ko 5 k . k, : A n k+ko 


g 


0 

Bk A ^k+ko — 5k A 5k+k 0 5 k +k 0 A Hk+ko 51 k +k 0 A W A Ek+k 0 / 




V 0 / 


where couples of self-energy parts are introduced E k = E(w, k) and E k = E(w, k), II k = n(w, k) and II k = 11* (a;, k), 
5 k = 5(w, k) and 5k = S(ui, k), Vk = V(oj,k) and V k = Ft(u;,k), and we use relations f? k = —Bk = — <B k+ko and 
5 k = — £ k = £ k +k 0 following from Eqs. (fTUl) . 

The general solution of Eq. m is quite cumbersome and we do not present it here. Green’s functions derived from 
Eq. (1121) in the spin-wave approximation (i.e., at zero self-energy parts) are presented in Appendix [A] (see Eqs. (|AljO . 
Their denominator has the form 


D (0) (w, k) = (w 2 - (e+ k ) 2 ) (w 2 - (e 0k ) 2 ) , (13) 

where 

(4) 2 = \ { E l + E l+ko -Bl- Bi +U A 2 Bl - 25 2 ) ± (14) 

= \ {El A A 2 +ko -Bl- Bl +ko A 2 Bl - 2 5 2 ) 2 

((-^k + -^k)(-£'k+ko -^k+ko) (^k — ^k) ) ((-^k -^k)(-£'k+ko -^k+ko) A(5 k A5 k ) 2 ) (15) 

and e± give energies of two magnon branches in the spin-wave approximation (i.e., the classical rnagnon spectrum). 
It is seen that dipolar forces split the spectrum into two branches as it was pointed out before . 19 ’ 20 It can be shown 
using Eqs. (HU, mi, and m that £p k are invariant under replacement of k by k + ko. 

To find magnon spectrum in the first order in 1/5 (that is denoted below as e^), one has to use Green’s functions 
m for diagrams calculation and to consider the first 1/5 corrections to the Green’s functions denominator that has 
the form up to a factor 


D (1) (w,k) = (w 2 - (e+ k ) 2 ) (u 2 - (e 0k ) 2 ) +D(w,k), (16) 

where f2(w, k) is a function linear in self-energy parts. The explicit expression for fl(u>, k) is given in Appendix [A] (see 
Eqs. (I ATI) and (I A5I) 1. It is also shown in Appendix lAl that fl(cu, k) is invariant under replacement of k by k + ko. As 
a result ef k (like £p k ) have the same form in the vicinity of k = 0 and k = ko. That is why we discuss below the 
spectrum only in the neighborhood of the point k = 0 (i.e., at k <C 1) bearing in mind that it has the same form near 

k = k 0 . 


C. Magnon spectrum 


Classical magnon spectrum given by Eq. HI becomes simpler in the limiting case of k <C 1. As it is seen from 
Eqs. ([TOP , one has to use properties of dipolar tensor at k <C 1 and k ~ ko to derive it. The dipolar tensor has the 
well-known form near the point k = 0 

dafU k a kft 

~3 P ~~) 


Q ap 

k = w 0 


k < 1. 


(17) 
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FIG. 2: (Color online.) The splitting of the magnon spectrum into two branches obtained in the spin-wave approximation (e^ k 
given by Eq. gsd) and in the first order in 1/S (e/ k given by Eq. (129m . Curves are drawn for S = 1/2, J = 1, and uo = 0.4. 
Inset shows a sketch of the sum of two dynamical structure factors iS^w, k) + S yy (ui, k). 


We obtain numerically using the dipolar sums computation technique (see, e.g., Ref.— and references therein) at 

k ~ k 0 

(/? k ^_ ko — Wo (pxx ^ ) ^ct/3 d - Cxykakfi (1 ^a/?)) , ^ 

C-xx » 0.051, (18) 

C X y ~ —0.055. 

In particular, it is seen from Eq. (fl8l) that Q k ^=0. It can be shown that higher order terms in powers of k in Eqs. m 
and (11811 do not contribute to the results obtained in the present paper in the considered orders in loq/ J and k. 

One obtains from Eqs. ED, ED, (ED, ED, and (fTSl) for the classical spectrum at k -C 1 in the leading order in 
Wo / J 


e 


± 

Ok 


Dkyjl - 2/c 2 L 2 (0 k ,< / ? k ) ± -^ii(^k)sin 2 9 k « Dk ^1 - k 2 L 2 (9 k ,ip k ) ± -^Li(</? k ) sin 2 6y 
Ti(v?k) = yj (1 + 12c xy ) 2 + 24(3 c xx - c xy ){ 1 + 18c ra + 6c xy ) cos 2 2<p k « sj 4.063 + 3.95 cos4<p k , 
L 2 (9 k , ip k ) = fcos 4 d k + cos 2 0 k sin 2 9 k + ^ sin 4 d k (7 + cos4^ k )) ss L 2 = 0.069, 


(19) 

( 20 ) 
( 21 ) 


where D = S^/2JJo = SJy/l2 is the magnon velocity, angles 9 k and ip k are taken in the spherical coordinate system 
with 2 -axis directed along the staggered magnetization and constants c xx and c xy are defined in Eq. (fl8l) . It should be 
noted that the classical spectrum lHUl) is gapless. The function L 2 (9 k , <p k ) is very smooth: its values lie in the interval 
[fl’ 12 ]- That is why L 2 (d k , ip k ) can be averaged over the angles and replaced by the constant L 2 for simplicity (see 
Eq. (12ip ) as it was done in Ref—. The spectrum splitting depends on the function Li(<p k ) which has the following 
properties: 


0.34 « 1 + 12 c xy = T 1 ( 77 / 4 ) < Li(ip k ) < Li(0) = 1 + 36c xx « 2.83 (22) 

and Li(<p k ) ~ 2y / 2| cos2i^ k |. Spectrum (fl9l) is plotted in Fig. [2] for a particular set of parameters. 

It should be noted that Eq. ED differs from the classical spectrum obtained in Ref.—. The origin of this discrepancy 
is in the fact that dipolar tensor components were found in Ref. 19 with the precision O(fc), whereas some quadratic 
in k terms contribute to £g k . In our notation, these are terms taken into account in Eq. (1181) (quadratic in k terms 
which are omitted in Eq. ED do not contribute to Eq. ED)- 

Magnon spectrum can be extracted from the dynamical structure factor (DSF) that is measured in neutron scattering 
experiment. In the spin-wave approximation, transverse DSF is determined by a linear combination of Green’s 
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functions m■ In particular, DSF S xx (ui, k) has the form at k ~ ko 

<S xx (w,k) oc Im(G(u;, k) + G(w, k) + F(u>, k) + F t (w, k)) 

= ^77 (<*(“> - 4) + s ( u + 4)) ( K l ~( 1 + 36 Cxx ^ C ° l ^^ 

+ ~ e 0 k) + <4 + e 0 k)) + (1 + 36c ra ) > 

where we use Eqs. (ED for Green’s functions in the spin-wave approximation. The sum «S xa: (u>, k) + <S w (u;,k) has 
the simpler dependence on angles: 

a j 

<S xx (w,k) +<S TO (w,k) oc 2ir —- (8(u> - e^ k ) + S(uj + ep k ) + 6(u - £g k ) + 6(u + e^)) , 

UJ 

where k ^ ko. It is sketched in the inset of Fig. [2j where we take into account that delta-peaks are replaced by 
Lorentzian functions due to the magnon damping derived below. 


III. THE GROUND STATE ENERGY RENORMALIZATION 


The classical ground state of the model 0 is continuously degenerate: the staggered magnetization can have 
arbitrary direction as it is seen from Eq. ©. It is well known that quantum fluctuations can give anisotropic 
corrections to the ground state energy selecting a limiting number of states (“order-by-disorder” effect). These 
quantum corrections are proportional in our case to sums over momenta containing components of the dipolar tensor 
Q J )' 3 and depend consequently on the direction of the quantized axis relative to the lattice. Thus, one should bear in 
mind in the subsequent calculations what is the easy direction of magnetization in the ground state. Using Eq. © 
for the biquadratic part of the Hamiltonian and Eqs. ED for Green’s functions, we obtain after tedious calculation 
the following anisotropic part of the first 1/S correction to the ground state energy E gs : 


A E gs 
N 


C = 


C, ,2 

JUJq , 2 2 , 2 2 , 2 2 \ 

C — \lxly + 1x1 z + Izly) 

J 1 

1EIa 

u q 


(Jo-J q ) 2 ((Q xx -Q^) 2 -4(0 x ^) 2 ) 


(Jo 2 - Jl) 


3/2 


0 . 0022 , 


(23) 

(24) 


where 7 j are direction cosines of the staggered magnetization relative to axes which are parallel to cube edges. 
Components of the dipolar tensor in Eq. (l24l) are taken relative to these axes. The constant G has been calculated 
numerically using the procedure of dipolar sums computation^. This computational technique is required because 
momenta q ~ 1 give the main contribution to the sum in Eq. EH) and one cannot use Eqs. (1171) and (fT51) . As G > 0, 
cube edges are easy directions for the staggered magnetization. 

It is shown in the next section that similar to FMs with dipolar interaction ^ ?-’ 1 - the fluctuation-induced anisotropy 
(1231) is naturally accompanied by the fluctuation-induced gap in the spectrum. Then, both the anisotropy and the 
gap have the “order-by-disorder” origin. 


IV. RENORMALIZATION OF THE REAL PART OF THE SPECTRUM 


Let us discuss renormalization of the real part of the spectrum stemming from diagrams of the first order in 1/S 
shown in Fig. |D Lines in these diagrams stand for bare Green’s functions introduced in Eqs. (fill) (see Eqs. (IA1I) for 
their explicit form in the spin-wave approximation). Each self-energy part arising in the Dyson equation (1121) receives 
its own contribution from the diagrams. 

As can be seen from results below, it is more convenient to discuss renormalization of the real part of the spectrum 
square for which we have from Eq. © 




Ref l (w = £p k + iS, k) 
2 v / dk 


(25) 


where the last term is given by Eq. (IA6I) . The Hartree-Fock diagram presented in Fig. ©a) originates from four- 
magnon terms © in the Hamiltonian. After simple calculations we obtain in the leading orders in k and wq/J for 
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the contribution to fl(u;,k) 


=F 


» (4) (4^) 

2i/rfk 




■Jo \ Jq Jq 


SN ^ Jo 

q 


(1 + 2A/”q) 


(26) 


-)xy 


J 0 S ^ (JO ~ Jq) 2 ((Qq - Q V q V ? + 2 (Q^f) „ , J 0 S ^ Qq + ko 

47V ^ / T 2 721^/2 1 + 7 qj + 27V ^ 


(^o 2 - -/ 2 ) 


7 2 — T 2 

J 0 J q 


(1 + 27V q ), 


where TVq = (exp(e q /T) — 1) 1 and e q = S^JJq — J 2 . One concludes comparing Eqs. (fl9l) and (1251) that the first 
term in Eq. (12611 leads to the well known renormalization of the magnon velocity D 


D^D 1 


2 57V 


Jo \l Jo Jq 

E— t ~- - (i + 2Ag 


Jo 


°hrsh^-9(i 


(27) 


where (/(x) is the Riemann zeta function and we assume that 5 ~ 1 (so that T -C D at T <C T)v). The second and 
the third terms in Eq. (126[) contribute to the spin-wave gap. 

The loop diagram shown in Fig- |TJb) comes from three-magnon terms ([5]) in the Hamiltonian. As a result of simple 
but tedious calculations we obtain for the contribution to the real part of Q(u>, k) in the leading orders in k and u>o/J 


=F Re 


^ (3) (4^) 

2\/3k 


3J 0 S 

2N 


E 


(J0-Jq) 2 (Q X q V ) 2 

(j 2 0~J 2 qf 2 


(1 + 27V q ) 



q 


0 

\]J 2 0 ~ J\ 


(1 + 27V q ), 


(28) 


where we set k = 0 under sums because the summation over q ~ 1 gives the main contribution. 

One obtains in the first order in 1/5 from Eqs. (1T51) . (1251) - (1251) and (1251) the following expression for the spectrum 
at 5 ~ 1 and k <C 1: 


c lk 


1 - 2 L 2 k 2 ± ^y-Ti(<p k ) sin 2 6 » k ) + A 2 , 


= \(Dkf 


where we imply the small renormalization of the magnon velocity m, 


(29) 


A = \/24 CSuq (30) 

is the gap in the spectrum, and the constant C is given by Eq. (l24l) . Notice that thermal corrections to the gap are 
negligibly small at 5 ~ 1. Eq. (1251) is plotted in Fig. [2] for a specific set of parameters. Spectrum (1251) has the following 
form in the two limiting cases: 


“lk 



1 - L 2 k 2 + 


A 2 


± 


w 0 


2 (Dk) 2 8 Jo 


Li(tpu) sin 2 6 > k 


k < A/D , 
k > A/D. 


(31) 


As it was done in Refs .- 13,17 for 2D and 3D FMs with dipolar forces, it can be shown that coincidence is not accidental 
of the numerical constants C in expressions for the anisotropic correction to the ground state energy (1251) and to the gap 

(1551) . Namely, the anisotropy in the Hamiltonian of the type Y2i {(^i) 2 (^i) 2 + (5f) 2 (5f) 2 + (5f) 2 (5f) 2 ) /S' 4 
(cf. Eq. (1251) '). where C is a positive constant, leads to the gap in the classical spectrum of the form (1551) if 5 1. 

As is seen from Eqs. m and m, the spectrum renormalization is very small at k 7 ?> A/D whereas quantum 
fluctuations change it drastically at k <C A/D. One has to take into account this renormalization when discussing 
the spin-wave damping. Then, we carry out below self-consistent calculations of the damping. Notice that such 
self-consistent consideration leads to the same result (1251) for the real part of the spectrum. 


V. MAGNON DAMPING 

It is well known that the magnon damping arises in Heisenberg non-frustrated AFs at T ^ 0 in the second order in 
1/5 and there is no damping at T = 0*2 Dipolar forces give rise to the finite damping at T > 0 in the first order in 
1/5 due to the three-magnon interaction (|5J) that leads to the loop diagram shown in Fig. Q2b). 























b) 


FIG. 3: (Color online.) Momenta q, k, and q— k (or k —q) of three magnons in the decay (a) and confluence (b) processes which 
are given by Eqs. 151 and 151, respectively. Components q|| and qx of q are also shown which are parallel and perpendicular 
to k, correspondingly. 


Contributions from the diagram presented in Fig. |ljb) to the imaginary part of each self-energy part contain delta- 
functions describing the magnon decay and the confluence of two magnons (it is clear from the explicit form of the 
bare Green’s functions given by Eqs. (ED). For a magnon with momentum k, one has 2 3 = 8 possible decay processes 
of the type 


which arise at any T and 8 confluence processes 


_ e ± _ e ± 
fc q fc k- 


= 0 


et - g 


fc k-q 


(32) 

(33) 


which exist at T / 0 only. Their contributions to the imaginary part of f2(u;, k) determining the damping are not 
equal and depend on T and k. Figs. Ha) andHb) illustrate Eqs. (l32l) and (1551) . respectively. The reader is referred 
to Appendix IQ for a detailed analysis of Eqs. (151 and (1531) . 

One obtains from Eq. (m for the magnon damping in the first order in 1/S 


r ± _ , Imn(h) = 4k + iS . k ) fo A \ 

k 4^4 • ( j 

As it was mentioned above, one has to carry out self-consistent calculations to find T^ due to the considerable 
renormalization of the real part of the spectrum by fluctuations at k < A. The general expression for rj^ and 
corresponding calculations are rather cumbersome. Fortunately, the results become quite compact in limiting cases 
which are important for our consideration and which we discuss below. 

We assume in this section that 5 ~ 1. As a consequence, it is implied that D ~ T'n and A ~ ujq. 


A. T = 0 


As it is discussed in Appendix [C] in more detail, the decay processes contribute at T = 0 to the damping of 
T”-magnon branch only whereas the “—’’-branch has infinite lifetime: 


Fk =0. 


As a result of simple but tedious calculations we obtain at k ^$> sin 0kand sin 0^ ^$> J 

Wg (1 + 12 C X y)~ 


F+ = h 

k J 2 276487rT 2 


xy> Ai(<^k) sin 4 dkcos 2 9k ~ fc^l.9 • 10 5 Li(ipi i ) sin 4 9^ cos 2 9^, 


(35) 


(36) 


and r,j( is negligible for other k and 6 k. Thus, one concludes from Eq. (1551) that r+ « e+ at T = 0. 

Temperature corrections to Eq. (1551) become important at T> Suj 0 sin 2 9^/k 2 . The order of magnitude of these 
corrections can be obtained by multiplying of Eq. (1551) by 


Tk' 2 


Su>o sin 0k 


■ In 


SiO o . 2 


sin 0k 


(37) 
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FIG. 4: (Color online.) Graphic of the function /(<pk) given by Eq. (14011 that determines dependence of the damping on 
azimuthal angle p k (see Eqs. (1551) . d39b . (l43l) . (l48l) . and (l49l l. 


B. T A 0 


Confluence processes give the main contribution to the damping when sin 0k Wf k ^ and sin 2 0k -C 

T^D: 


rj = 


wg T(l + 12^ a^— ^ f 

kJ D 2304t r V ; V 

/(+ k) 


/Ok) 


4(1 +18c xai + 6Cxy) 


sin 2 2#k, 


r k rk 4(1 + 18c ra + 6c xy ) - /(+k) ’ 


where the non-negative function 


/(+k) = 


(1 + SOCjEa; i/1 ( + k ) ) (^1 (+k) 1 12Ca;y) 


Li(tpk) 

is introduced (cf. Eq. ([22]) ) which graphic is shown in Fig. [4] and 
A = 

144 L 

288L 2 A 2 6912CL 2 


144L 2 J+ l( *’ k)8i1 ' 2 " k ” 0 ' 1 7t Llfek) ™ 2 9k ’ 


5 = 


1.05 


S 2 ujqLi(p\A) 2 sin 4 0k STi^k) 2 sin 4 0k S'Li(+k) 2 sin 4 0k 


(38) 

(39) 

(40) 

(41) 

(42) 


Eqs. (1551) and (1551) are valid for B < 1. Momenta of summation g > fc give the main contribution to Eqs. (1551) and 


For larger momenta, one obtains when k 3> sin 0k y/^f, sin 0k J and T Dk 


r+ = fc 2 ^i- (i + i 2c ^) 2 sin2 


J D 


768t r 


20 k 


x f ln ( s \fJwo k gin 0 

V V A 

r - _ ^ 2^0 T (1 + ^°xy) 2 ^ 


1 - 12 c x 


xy 


J D 


768tt 


1 + 12c ; 

Vb n 
l - VT^b 


i - 


/Ok) 


4(1 +18c 

XX ^C-xy) 


In 


Vb 


i - VT^b 


sin 2 2^k, 


(43) 

(44) 


where decay and confluence processes lead to the first and to the second terms in the last brackets in Eq. ©, 
respectively, and confluence processes determine . 

It is seen from Eqs. (1551) (1551) and (1551) (1551) that the damping T^ decreases as k 2 upon k decreasing down to 
k ~ VW w hich region this decreasing changes into rising that has the form 1/fc 4 . This rising takes place up to 
k ~ A/D near which point the increasing turns into a rapid fall due to the gap in the spectrum. Thus, T^ has a 
peak at k ~ A/D which height is of the order of uiqT/D and the ratio r^/e^ is proportional to T/D <C 1 at the peak 
position (see Fig. [5]). Thus, one concludes that magnons are well defined quasiparticles in quantum AF with dipolar 
forces at T <C T/v. 
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FIG. 5: Sketch of the ratio r^/e^ at T SyflToo. Asymptotic at k y/uo/J is also shown. The peak hight at k ~ A/D 
is proportional to T/D -C 1 for quantum spins and it is given by a constant for S > 1 at D < T < Tn- The peak height in 
quantum AF can reach a value of the order of unity if the gap is sufficiently reduced by the anisotropy competing with the 
dipolar one (1231) (see discussion in Sec. IVI SI . 


VI. POSSIBILITY OF THE MAGNON BREAKDOWN 

We obtain in the previous section that the damping rising upon k decreasing stops at k ~ A/D due to the gap 
in the spectrum. One infers that a reduction of the gap value could keep the damping increasing and lead to the 
long-wavelength magnon breakdown. We discuss in this section two possibilities of the gap decreasing. First, we 
consider large spins S 1. As the gap (15U1) obtained above is of the next order in l/S as compared to the bare 
spectrum, the gap value can be reduced by increasing S. The second way to decrease the gap value is to take into 
account a magnetocrystalline anisotropy competing with the dipolar one (1231) (i.e., the magnetocrystalline anisotropy 
favoring cube space diagonals rather than cube edges). In this case, a negative contribution arises to the gap (15(1 
which can decrease the gap considerably (see Appendix [D] for a more detail discussion of the possibility to reduce the 
gap value in this way). 


A. Large spins 


It is well known that TV ~ S 2 J whereas D ~ SJ. Then, the temperature can lie in the range Tn T D 
for large enough S. Besides, quantum corrections to observables decrease upon S increasing and they die out in the 
limit of S —> oo. In contrast, ratios of temperature corrections to the bare values of observables contain powers of 
T/S 2 J ~ T/Tn so that temperature corrections remain finite in the limit of classical spins (see, e.g., RefJ^ for detail 
discussion of this point). As a result, to calculate the gap at T N T D and S 1, one can replace 7V q by T/e q 
in Eqs. E51) and (E51) and discard T-independent terms. We obtain in this way for the gap instead of Eq. (1501) 

A|> = C»S 2 u; 0 2 ^, (45) 

where the constant C> is defined as (cf. Eq. EH) ') 


_ 3J 2 1 v (JO ~ Jq) 2 ((Q^ ~ Q V g V ) 2 ~ 4(Q q y ) 2 ) 
N ^ ( J o 2 - J 2 ) 2 


(46) 


The summation over q ~ 1 gives the main contribution in Eq. (1461) . 

Damping estimation leads to Eqs. (1551) and (1551) for sin$k\/^- » k and Tn T D, where now B = 0 

and A is a constant of the order of unity. These expressions give for ratios ~ WqT//c 2 J£) 2 . Then, one obtains 

using Eq. (1551) near the peak position at k ~ -gp and at fixed ipu and 9^ 

p± 

—j- ~ const. (47) 

e k 

Thus, we demonstrate a breakdown of a small fraction of long-wavelength magnons (with k ~ and sin 20k ~ 1) 
for S 1 at small temperature T <C Tn- 
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B. Competing magnetocrystalline anisotropy 


We assume now that the gap A in the spectrum (l30l) is decreased by the competing magnetocrystalline anisotropy 
so that A -C VScoo- Counterparts of Eqs. (1551) and (1551) which are valid at S ^t sin0k » k sin 0k and 

D T ~q 2 — tt ■ '> a have the form 

S 2 uJok sin^ 9^ 


r£ 

r k 


(T\ 4 n 3 (l + 12c xy ) 2 / f(<p k ) 

kJ\D) 5760 ^ 4(1 + 18c xx + 6c xy ) 

~ jj 6.1 • 10- 4 (1 - 0.16/(v? k )) sin 2 20 k , 

r +_ fW k) _ 

k 4(1 + 18c 

XX + 6 C X y) - /(95k) ' 


sin 2 20k 


(48) 

(49) 


One concludes from these equations that the damping becomes of the order of the real part of the spectrum near the 
peak position at k ~ A/D if the gap value is decreased so that the inequality A < uj 0 (T/D) 2 -C ujq holds. 

We expect that there is a small chance of success to find a cubic AF in which the magnetocrystalline anisotropy 
cancels almost completely the dipolar gap. For instance, the anisotropy in the most perfectly isotropic cubic AFs 
TlMnF 3 (Ref/ 22 ) and RbMnF 3 (Refs^r— ) competes with the dipolar one being of the same order of magnitude. But 
it turns out to be slightly greater in these substances than the dipolar anisotropy so that the easy directions are space 
diagonals of the cube. The resultant gaps in these materials turn out to be even slightly greater than the dipolar gap 
given by Eq. (1501) . 

However a way was proposed to change gradually values of the anisotropy (and the gap) in TlMnF 3 (Refj^) and 
RbMnF 3 (Ref. 2 -) by replacing a tiny amount of Mn 2+ ions by Co 2+ . As Mn 2+ ions are in spherically symmetric 
states with L = 0 and S = 5/2 in these compounds, the magnetocrystalline anisotropy is tiny so that the anisotropy 
field favoring (111) directions is equal to several oersteds only. In contrast, Co 2+ ions have L ^ 0 and S = 3/2. 
As a consequence, the effect of spin-orbit interaction is much more pronounced: the anisotropy field selecting (100) 
direction is four orders of magnitude larger than that of Mn 2+ . Thus, two single-ion anisotropies on Mn 2+ and Co 2+ 
ions compete in mixed compounds TlMni_ 2 .Co a; F 3 and RbMni_ x Co x F 3 . Due to the great difference between the 
anisotropies magnitudes on Mn 2+ and Co 2+ , a very small x = a;* is required to change the easy direction of the whole 
sample from (111) to (100): x» ~ 0.0004 and a;* ~ 0.00034 for TlMni_ x Co x F 3 and RbMni_ x Co x F 3 , respectively. 
The gap value is reduced considerably at x ~ cc*. 

It has been shown recently^ that states near the rnagnon band edges can become localized in disordered systems 
with gapped spectrum. However, according to estimations made in Appendix [Dj states with k x(A!/ J) 2 ~ 
10~ 6 A -1 , where A! is the value of the magnetocrystalline anisotropy on Co 2+ , remain propagating in materials under 
consideration. Besides, the rnagnon damping due to the scattering on impurities is negligible at such k. On the other 
hand, Eqs. (1551) and (1551) predict the rnagnon breakdown at k < 10 -3 4-10 -4 A -1 due to the rnagnon interaction with 
each other if the gap is reduced considerably. Notice also that momenta of summation q < 10~ 6 A -1 are inessential in 
the calculations leading to Eqs. (1551) and (1551) . Then, TlMni_ x Co x F 3 and RbMni_ x Co x F 3 can be suitable for testing 
of our predictions. 

Unfortunately, it would be difficult to carry out corresponding experiments because the characteristic values of 
momenta of overdamped magnons are quite small being of the order of 10 -3 -f-10 -4 A -1 . However, bearing in mind 
recent progress in neutron spin-echo techniqu o 15 ’ 16 , we hope that the corresponding measurements will become feasible 
in the near future. It is also possible that more suitable substances will be found which have larger values of momenta 
at which the discussed anomalies arise in the damping. 

It should be noted also that our conclusion about suitability of the mixed compounds for the observation of the 
rnagnon breakdown is based on estimations made in Appendix [D] in the first order in x. As soon as defects change 
the bare spectrum considerably at x ~ x* and fcCl, these estimations must be used with caution. In particular, one 
cannot fully exclude the possibility of great spectrum change by terms of higher orders in x. It is difficult to analyze 
the whole series in x but we point out that all the expected contributions are small as compared to those of the first 
order in x due to the smallness of all kinds of anisotropy in comparison with the exchange constant. 


VII. CONCLUSION 

To conclude, we discuss rnagnon damping in Heisenberg AF on a simple cubic lattice with dipolar forces at small 
temperature T <C Tjy. In accordance with previous results, it is demonstrated that dipolar forces split the rnagnon 
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spectrum into two branches. The classical gapless spectra of long-wavelength magnons into two branches are given 
by Eq. (1191) . It is found that quantum and thermal fluctuations modify the spectrum considerably near points k = 0 
and k = ko : the gap A (see Eq. (151)1) ) appears in the spectrum. The gap is accompanied by anisotropic corrections to 
the ground state energy (ESI) which make cube edges easy directions for the staggered magnetization. These effects 
are of “order-by-disorder” nature. The renormalized spectrum of long-wavelength magnons is given by Eq. m- 
It is shown that magnons are well defined quasiparticles for all k at S ~ 1 and the ratio T k /e k has the peak 
at k ~ A/D which height is proportional to T/D -C 1 (see Fig. [5J. We discuss some possibilities of observing a 
phenomenon contradicting expectation of the quasiparticle concept: the breakdown of some part of long-wavelength 
magnons. In particular, it is shown that T k /e k ~ const (at fixed ipk and 9 k) near the peak position when S 1 and 
T n / S -C T -C T n . It is also shown that a single-ion anisotropy which competes with the dipolar one (E51) reduces 
the gap value enhancing the peak height. The peak height can reach a value of the order of unity for sufficiently 
small gap that signifies the breakdown of long-wavelength magnons with momenta lying near the peak position. The 
gap can be decreased and the magnon breakdown can be stimulated also by replacing of a small amount of magnetic 
atoms by those with single-ion anisotropy competing with the dipolar one. We argue that this effect can be observed 
in TlMm_ x Co x F 3 and RbMni_ x Co x F 3 at x ~ 0.0004. 


Acknowledgments 

This work is supported by Russian Scientific Fund Grant No. 14-22-00281. 


Appendix A: Green’s functions and general expression for the spectrum 


Solution of Eq. m has the following form in the spin-wave approximation (i.e., with zero self-energy parts): 
(iO + El/) (w 2 + -Bk_|_ko — 4c+k 0 ) 2Sk<?k^k+k 0 + <?k ( w — ^k+k 0 ) — ^k ( w A ^k+k 0 ) 


G(w,k) = 

^(w.k) = 


-£k (t o 2 + B*-E? 


(c ^ 2 - (e+k) 2 ) (^ 2 - (e 0k ) 2 ) 
k+k 0 — -'Ac+ko) — 2Bk£kEk+k 0 + Bk+k 0 ( B k + £ } 


( w2 ^ (4) 2 ) ( w2 - ( e 0k) 

+ , . £k (Sk+ko ( w + Ek) + -Bk (to — E k +k 0 )) ~ Bk (Sk^k+ko + (w + Eh) (to — Eh+ ko) + ^k ~ ®k) 

F' (to, k) = - 7 -- 7 T-, 

( w2 - (4k)) ( w2 - ( e ok)) 

Bk (Bk ( to + -Ek+ko) — Bk+k 0 ( w + Ek)) + £k ( — BkBk+ko + (w + Ek) (to + Ek+ko) + £k ~ B>k) 


(Al) 


G(w, k) = 


( w2 - (4k) 2 ) ( w2 - ( e ok) 2 ) 


where energies ej k of the two magnon branches have the form (flTl) . By setting ioq 
to the Green’s functions of Heisenberg AF (see, e.g., Ref^): ^(Wjk) = Q(to, k) 
F^(w,k) = Fj(w,k), where 


0, one leads from Eqs. ED 
0, G(w,k) = G c (to, k), and 


G c (to, k) 
Ft (to, k) 


S Jq uj 

u>*-s*(j$-4y 

SJk 

oj*-s 2 (j 2 -Jiy 


(A2) 

(A3) 


The explicit expression for f2(w,k) introduced in Eq. (flfil) has the form 

fl(w, k) = S£l(to, k) + 8£l(to, k 4- ko) 


(A4) 








13 


where 


k) = (B k ( w 2 + E 2 +ko - El +k0 ) + 25k5kE k+ko - E k+ko (Bl + 5 2 )) (n k + 4) 

- (E k (u, 2 + S 2 +k0 - E 2 +ko ) + 25 k 5 k S k+ko - E k+ko (5 2 + 5 2 )) (S k + E k ) 

+ w ( w “ + -®k+k 0 ~ ^k+k 0 + &k — ®k) (^k — S k ) 

+ W (5 k (Ek — E k + ko ) + 5k (Ek + E k + ko )) (j^k + E k ^ 

- uj (£> k (E k — B k + ko ) + 5 k (-Ek + E k + ko )) (5 k + 5 k ) 

+ (5 k (Bk+ ko E k — E k E k+ko ) + £> k (cj 2 + -B k S k+ko — E k E k+ko + 5 k — B k )) (v k — V^j 
+ (Bk (E k E k+ko — -B k+ko E k ) + 5 k (uj 2 — B k B k+ko + E k Ek+ ko + 5 k — B k )) (5 k — 5 k ) . (A5) 


One obtains for the first 1/5 corrections to f2 (uj = 6g k + id, k) at k -C 1 in the leading order in k and ujq/J 


T 


f 1(uj = e. 


o k 


id, k) 


2v / d k 


SJp 
2 

SJ 0 


+ 


lT(l + 36c ra ) ii ( (( J k k ) ) ) ( 

1 T (1 + 36c xx ) 


'i(<£ k ) J 

cos(2 <p k ) 


n k + if k — s k — 


s k ) 


Li(iph 


(n k+ko 


+ n k+ko + E k + ko + E k + k 


± iSJo- 


1 T 12ca;y sin(2^j k ) 

2 Li(ip k ) 


(E k 


-vl+ S k - Sk + B k+ko - 'P k+ko 


- k 


SJp 

2>/3 


(E k - S k + E 


k + k Q E k _|_ k() 


3 ) 


uc t 1 + 36c xx cos(2<y9 k ) ^ ^ ^ 

^/COJq ^ /7T r ( . _ \ V^k -^k ^k+ko ' ^k+ko ) 


=h ikSJo 


2-\/3 Ei(^k) 

1 + 12c gy sin(2y; k ) 
2 V / 3 El(<y3 k ) 


(5 k + Sk — 5 k + ko — 5 k + ko ) 


+ k 

+ k' 


2 SJo 
~24~ 
2 B Jo 

~24~ 


1 ± (1 + 36 c xx ) 


C0s(2<^i< 


=F ik 2 SJ 0 


— 1 ± (1 + 36c xx ) 

1 + 12 c xy sin(2</? k ) 
24 El(v?k) 


Li(<p k ) 

cos(2t^ k ) 


)(' 


n k + n k + E k + s k 


Li(ipk) 


(n k+ko 


n 


k+ ko 


^k+lco E k -f kc 


(A6a) 
(A6b) 

+ 5 k + ko — 5 k + ko j|A6c) 
(A6d) 
(A6e) 
(A6f) 
(A6g) 
(A6h) 


(' Pk 


- 4 - 5 k +5 k +E k+ko - 4 +ko 


5 k+ko + 5k+^^6i) 


where self-energy parts are taken at ui = e^ k + id. All terms in Eq. (IA6I) contribute to results presented in the main text 
for the damping due to decay processes, while only terms (IA6l) fal-fcl give the leading contributions to the damping 
due to confluence processes. Expressions for combinations of self-energy parts which arise in Eq. CM} are presented 
in Appendix [Bl 


Appendix B: Expressions for self-energy parts 

In this appendix, we present expressions for some combinations of self-energy parts which arise in Eq. (IA6I) . Only 
contributions are shown below which are of the first order in 1/5 and which originate from the loop diagram depicted 
in Fig. [T}b). To make all expressions more compact, we move arguments of Green’s functions to subscripts and 












introduce the following notation: k = (w, k), q = (u) q , q), and ko = (0, ko). 


n k 


- s k - x k = 


2 ({Fk-q 'Ff.-g + Gk-q Gk-q){Fk 0 +q Fk 0 +q Gk 0 +q + 0k 0 +q)Q]i-qQq 

q 


w„,q 

+ (-Pfc+fco-g — Fk+ko-q Gk+ko-q — Gk+k 0 -q)(Fq ~ F^ — G q + G 9 )Q k + ko _ q Q 

■^X2:2 


>212: 

q 


9 )^ +q Qr +k0 


- Fl +ko _ q - G k+ k 0 - q - G fc+fco _,)g**g£ q 

^yz 


+ {Fk+k 0 -q + -Ffc+fc o-q ~ Gk+ko-q - Gk+k 0 -q){F q + Fj + G q + G q )Q J Cj 
F{F k -q + F\_ q Gk-q — Gk-q){Fk 0 +q + F kq+q + Gk 0 +q + Gk 0 +q)Qc l Q ko +q 

-2 i(Fk-q + Gk-q)(Fl o+q - Gko+qW^Qr+Uq ^k+ko-q + Gk+ko-qM J - ^g)Qk+ ko - q Qk+ ko 

+2i(Ft +G,)(7- fc _, -a fc - g )O^Qr +ko - 2 i(Fl +k0 _ q - G fe+fco _ 9 )(F feo+? + G ko+qWJ QV+u 
+HFl 0+q +Gk 0+ q)(F k+ ko-q Gk + ko-q)Ql S Q^ +ko ~ ^lq ~ Gk-q){Fq +G q )Q™ 

+ 2 i{Tk+k 0 -q — Gk+k 0 -q){Fq + Gq)Q q Q k + k() + 2 l{Fk-q — Gk-q){F ko+q + Gk 0 +q ) 
~^{Fk+ko-qF > l 0 J rq + Fk—qF q + Gk+kq-qGk 0 +q + Gfc_ 9 G 9 )Q k+ko 
Fk 0 +q + F kq+q + Gk 0 +q ~ Gk 0 +q)(Fk+k 0 -q 
+ 2 {—Fk+ko-q + S/fc+fco —<3)(-^fco-i-q ~ ^fco+«)Q k + ko Q k - q 
+ i(Tk-q + P k _q - Gk-q Gk-q)(~F q + F* + Gq - GJQ-Q^ko-q 

+2(-F fc _ g + G k - q ){Fl - G 9 )gr +k0 Qr +k0 - q + 2 (-Fk+ko-q + Gk+k 0 - q )(J r l +q Gk 0+q )Q y k z +ko Q« z 

Fi{—Fk 0 +q + F ko+q + Gk 0 +q — Gk 0 +q)(Fk+k 0 -q + F k+kq _ q — Gk+k 0 -q ~ G k+k 0 -q)Qq Qq 
+i(Fk 0+ q + -^fco+q — Gko+q ~ Gk 0 +q)(~Fk-q + F k _ q — G k -q + Gfc-q )Q k -qQ ko + q 
+*( — • Fk+ko-q + F\ +ko _ q — Gk+k 0 -q + Gk+k 0 -q){Fq + — Gq — Gg)Q k+ko _ q (2 kc|+q 

+i(F k -q + J-L g - Gk-q - Gk-q)(-Fq + i* + Gq - G g )Q*^ +q 
+*(-?~fc+fco-q + F\ +ko _ q — Gk+ko-q ~ Gk+ko~q){~ Fq + Fj — G q + G g )Q q Q ko+q 
+*( — Fk 0 +q + -^fc 0 + 9 — Sfco+9 + Gk 0 +q)(Fk-q + F^_ q — Gk-q — G fc-<7)Q ko + q Q ko + q 

+ 2(F^ o+? — Sfco+gX - Fk+ko-q + Sfc+fc 0 -9)(3 k+ko Q ko+q 

+ 2(-^|_j_ feo _ g — Gk+k 0 -q){ — Fk 0 +q + Sfc 0 +«) ( 3 k + k o < 3 k o+q ~ Gk-q)( — F q + G q )Q k+ko Q ko + q 

+2(Ft - G,)(-F fe _ ? + G fe - 9 )Qr +k0 ^: +q + 2(-F fe _ 9 + G fe _ ? )(Ft - G 9 )gr +ko Q^ +q 

— (-^fc+fco-? + F k+ko _ q — C/fc+fco — <3 — ^fe+feo-9)(-^'fco+g F kq+q — Gk 0 +q ~ ^fco+9)Q k - q Q ko + q 
-{F k -q + F k _ q Gk-q ~ Gk-q){F q + F^ q - G q - G 

— {Fk+ko-q + F k+ko _ q — Gk+ko-q ~ Gk+k 0 -q){Fk 0 +q + F kQ+q — 

— {F k - q + F k _ q - Gk- q - Gk-q){F q + F^ q -G q - G ? )g^f +q ) . 


' F Fl o +q Gk 0 +q Gko + q)Ql g k0 +q 
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Ilk + l4 “I ^k + Sk — 

4“ Fk+ k Q — qF q “h Gk—qG kq+q "I Gk+kq— gG g )Q k 

w,,q 

+2(J- fc _ ? + <? fc _ 9 )(J-| 0+g +^ 0+9 )Qk a Qk- q + 2(F fe+fc0 _, + Gk+k 0 -q)(Fq + G,)QrQk Z + ko- q 
+2(J-I 0+9 + e fco+ ,)(7] fc _, +a fc - 9 )QrQq Z + 2(4- 9 + Gk- q )^k 0+q + Gk 0+q )Ql z Ql z 
+2 (Ft + Gq)(Fk+k 0 -q +Gk+k 0 -q)Q]< Z Qc l Z + ^4+kq-q + G k +k 0 -q)(F q + G q )Q^ Z Q^ 
+2(F k+k0 - q + G k+ko . q ){Fl + G g )Qk Z Qq Z + 2(7- fc _, + ^- 9 )(4 0+9 + ^*o+,)QrQ^ +q 

+ (-Ffc-g + -4c- g + + <?fc-g)(-Ffc 0 +<7 + Fk 0 +q + ^o+g + 5feo+g) ( 3k-qQq~ 

+ (-Ffe+fco-g + F k+ko -q + G k+ko - q + G k +k 0 -q)(F q + + Gq + G' (? )Qk+k 0 -q ( 3q' 

+ (-Pfc+fco-g + Fk+ko-q + G k+ k 0 -q + Gk+k 0 -q)(F q + F £ + G g + G g )Qq“ 2 

+ -4c_ 9 + 0fc-g + Gk-q){Fkq+q + -^fe 0 +g + 5fc 0 +g + Sfc 0 +g)Q q Qk 0 +q 
+2i(— Fk+kq-q + G k +k 0 -q){F ko+q + {?fc 0 +g)Qk Qk-q + 2i[—F k - q + Gk-q)(Fq + G g )Qk Qk + k 0 -q 
+*(-7 r fc 0 +g + F ko+q + Gk 0 +q + Gk 0 +q)(—F' k + ko - q + -F^+fc 0 _g + Gfc_ \-k 0 - q — G k J rko - q )Q ^ Q k _ q 
+i(~F k -q + T k _ q + £fc-g - G k -q)(F q + + G q + G q )Q^ Qk+ ko -q 

Fi( — F ko +q + F ko+q + Gkq+q — Gk 0 +q)(F k+ko _ q + F k _^ ko -q + G k -\- k q-q + G/ c +fc 0 _ g )Q q Qq 

+2i(Ft - G,)(J- fe+fco _ g +Gk + k 0 - q )Ql z Q V C+^ + 2 *(4-, + G fc _,)(-J- fc0+g FGk 0+q )Ql z Q y C+^ 

+2i(F ko+q — Gko+q)(F k -q + Gfc_ g )Qk Qk 0 +q + 2 i[F k -q + Gk-q)(F q — G q )Q k Qk 0+q 
+2i(F k+ko — q + Gk+k 0 -q){~Fq + G g )Qk Qko+q + 2*C F k +k 0 -q + G k +k 0 -q) (-^fc 0 + g — Gk 0 +q)Q]i Qq 
Fi(—F ko +q + F ko+q — Gkq+q + Gk 0 +q){F k -q + -Ffc_ g + Gk-q + Gfc_ g )Qk- q Qk 0+q 
+*(-Ffe+/c 0 -g + -F"fe-|-fc 0 _g + Gk+ko-q + 5fc+feo-g)( — F q + F^ — G q + G g )<5 k+ko _ q (5^ o+q 
+i(J- fc _ g + + S fc - g + Gk-q){-Fq + F\ + G q - G g )Q^ +q 

+*(-Ffe+/c 0 -g + ■F'fe-(-fc 0 _ g + Gk+kq — q + Gk+kq-q)(~F q + F^ — G q + G g )Q q Q ko+q 
+*( — -^fco+g + Fk 0 +q ~ Sko+q + Gk 0 +q){F k -q + + G k - q + Gfc_ g )(5 ko _)_qQ ko _|_ q 

— (-Ffc+feo-g — F k+kg _ q — G k +k 0 -q + Gk+k 0 -q)(Fk 0 +q ~ F kg+q + Gk 0 +q — 5fco+g)Qk-qQk 0 +q 
~(Fk-q - F fe t_ 9 - Gfc_ g + Gk-q)(Fq - F^ q +G q - S,)Q"* +ko _ q Q^ +q 

— {Fk+ko-q + F k+kg _ q + Gk+ko-q + Gk+kq-q) (Fk 0 +q + F kg+q — Gk 0 +q ~ Gk 0 +q)Qq Qk 0 +q 
~(Fk-q + Fl_ q + G k - q + Gk-q){F q + F+ - G, - G ? )Q^f +q ) . 


(B2) 
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~ F k — 5k + 5k — 

- 2^ T E (Vj + - Zk- g )Q&M z ~ 2 (Fl +ko _ q - G k+ko - q )(T ko+q + Qk 0+q )Ql z +ko Q 

iO q ,C{ 

+2(E 0+9 + G k0+q )(F k+k0 . q - G k+ko - q )Qll ko Q™ - 2 {Fl_ q - G k - q )(F q + G q )Q^ +ko Q^ 

+(F k . q - l\_ q + g k - q - g k - q ){F q + Ft + G q + G g )Qr- q Oq 2 

+ (-?fco+g + -F ko+q + g ko + q + g ko +q){F k + ko — q — F^ +ko _ q + G k+ko - q — Gk+k 0 -q)Q\t,+^ 0 -qQ. 

+ (^fc-g — -^l- 9 + g k - q ~ gk—q') (,F q + Ft + G g + Gq)Q, 


(B3) 


1012 : 

q 


-\xz 

*q 


-\xz2 

/q 


+ (-5feo+9 + -5fc 0 + 9 + g ko +q + 0fe o +(j)(-Pfc+feo-9 F^ +ko _ q + Gk+ko-q Gk+k 0 -q)Qq Qk 0 +q 
+4 i(F k + ko - q F ko+q + F k — q F q + gk+ ko - q g ko + q + Gfc-gG g )Qk+k 0 Qk+ko 
+2 i(F k - q + G k -q)(Fl + G 9 )Q£l q g^ +ko + 2*(F fc+fco _, + g k+ko - q )(H 0+q + Sk 0+q )Q^ 0 ^Q v k+ko 
+ 2 i(F k . q + Gk-q)(F} +G q )Ql Z Ql Z +k0 + 2 z(F k + k 0 -q + & +feo - g )(^ 0+g + Sfco + ^+qQ&ko 
-i(F k+ko - q + Fl +ko _ q - Gk+ ko -q ~ Gk+k 0 -q)(F q + fJ + G g + G 9 )g q ~g k l q 

+2(^ r fc+fe 0 - g - g k + ko - q )(F jt + G ? )Q k+ko Q k _ q + 2 (F k - q — Gk- q )(F ko + q + <?fc 0 + g )g k+ko g k+ko _ q 

—i{F k -q + ~ — g k - q ){F ko + q + F ko+q + g ko +q + Sfc 0 +<j)Q q Qk+k 0 -q 

-i(Fk-q - Fl_ q + <? fc _ g - g k - q )(F k0+q - F fe t o+g - g k0+q + g k0+q )Qi z Qi z 

—2(F k - q + G k - q )(Fl o+q — Sfco+qjQk+ko^q' — 2(Ffc + fe 0 _ g + g k+ko - q )(F i ’J — G g )g k+ko g ko+q 

F2i{F kq+q — g ko +q)(F k +ko-q ~ ^fe+fco-«)Qk+k 0 Qk 0 +q + 2 *(F k+fco _ ? — ^/fe+fc 0 -g)(-?fco+g — ^fco+?)Qk+k 0 Qk 0 +q 
+2*(Ft - G q )(F k - q G k -q)Ql z + k0 Ql z 0+cl + 2i(Fl_ q - G fc _ 9 )(F 9 - GJQZ^Q^ 

-i(F k - q - F\_ q + G k - q - Gk- q ){F q - Ft + G q - G 9 )0£l q Q^ +q 

i(Fk+ ko -q — F k+ko _ q + g k +k 0 -q — g k +k 0 - q )(Fk 0 +q — F kg+q + Sfc 0 +q — ^fco+<?)Qk+k 0 -qQk 0 +q 

-UK-* ~ Ft_ q + Gfc_ 9 - G k - q )(F q - Ft + G, - S,)Q“Q^ +q 
-~i(F k +k 0 -q - F k+kg _ q + Gk +k0 -q - G k +ko-q)(Fq Ft - G, + G g )Q-g^ +q 
i{Fk+k 0 -q ~ -Fk+ko-q "t" gk+ko-q — 5fe+fc 0 - 9 )(-Ffe 0 + 9 — F kg + q + 5fc 0 +(j — ^fco+<?)Qk 0 +qQk 0 + q 
+ (-?fe+feo-4 + E+fc 0 -g — g k +k 0 -q — g k +k 0 -q){~Fq + Fj — Gq + G q )Q k _ q Q ko+(i 
F{—Fk 0 +q + F ko + q ~ g k0 +q + g ko +q)(F k -q + F\_ q — Gk-q — G k -q)Q k + ko - cl Q k o+ cl 
F{Fk 0 +q + F’ko+g — g ko +q — gkq+q) ( F k —q + F k _ q — Gk-q + Gk-q)Qq g^g+q 
~k+k 0 -q — g k +k 0 -q + g k+ko-q){F q + Fj — G q — Gq)g ko+q ^ . 


+ (—Fk+ko-q + F t 
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’^Sk-Sk 


i T E(-« 


(B4) 


—Q 


<*>q, q 


&-,)(*? - G,)Q£‘Qk- q - 2{Fk+k 0 -q + G k+ko - q )(P ko+q - <5* 0+g )QrQkW 

+2(J- fc _, - &_,)(** + G 9 )gr^r + (-fife-, + 4- g + &- 9 + - Ft - G 9 + G g )Qr- q Or 

+ (“^fco+9 — ^kq+q ~ Gk 0 +q + Gk 0 +q){Fk+k 0 -q + -Pfc+fe 0 _ g 


+ l- / 'fco+9 — Sk 0 +q _ t/fc 0 +g "I" yk 0 +q){ r k+k 0 -q + -Ffc_|_fc 0 _ g + Gk+k 0 -q + G k +k 0 —q )Qk+k 0 -q^q 
+ (Fk-q — Fk-q ~ Gk-q + Gk-q)(F q + Fjj + G q + G q )Q q Z2 + 2(F k+ k 0 -q ~ G k+ k 0 -q)(F ko+q + 
F{Fk 0 +q + -Fk 0 _)_ 9 + 0fc o +9 F Gk 0 +q){Fk+k 0 -q ~ F k+kq _ q — Gk+k 0 —q + G k +k 0 -q)Qq Qk 0 +q 
+4*(Ffe_ 9 Fj 0+g + Fk+ka-qF^ + Gk-qGko+q + G k +k 0 -qG q )Q^~ 

+ HFt 0 + q + Gk 0+ q){Tk-q+Gk-q)Q^Ql Z + 2 i(P k _ q + G k-q) (F ko +q+ G k 0 +q)Q^ Q£* 

+ 2*(Ft + G q )(Fk+k 0 -q + Gk+k 0 -q)Qq Z QlG + 2*(F fc t +fco _, + G k+ k 0 -q ) (F, + G q )Q% 
+2i(F k+k0 -q - G fc+feo _ g )(Ft - 3,)Q£’Q"iq + 2i{F k + ko - q - G fc+fco _ 9 )(Fj - G 9 )Q£^ +q 
-i(F k +k 0 - q - Fl +ko _ q - Gfe+fc 0 _ g + G fe+fe0 _ g )(F, - Fj -G q + G q )Q™Q ^_ q 
+2z(F fe _, - ^- 9 )(^ 0+9 - ^ 0 + 9 )<3k Z <3k+ko-q + ^ &-«)(^£,+ g - G ka+q )Ql z Ql z 

, t , 7T \f-r t ~t i ~n \r\xzr\yz 


,xz 


Gk 0+ q)Ql Z Ql Z 0+ q 


—i(F k - q — F k _ q — G k —q + Gk-q)(Fk 0 +q ~ F\ a+q — Gk 0 +q + G ko +q)Qq i Qk+k 0 -q 

— i(F k ~ q — T k _ q _ ... , . . 

—i(F k - q + Fl_ q + G k - q + Gk~ q )(F q + Fj — G q — G q )Q^ 


^q+Gk 0+ q)Ql Z Q^ 


k-qQko+q 
rt 


iys 


n yz 

q^ko+q 


i(^k+fco- 9 -^k+ko-q F Gk+ko-q + G k+ko-q){F ko +q + F kq _y q Gk 0 +q Gk 0 +q)Q]i+k 0 
-i{Fk-q ~ F k _ q - Gk- q + Gk-q)(F q - Ft + G q - 5,)Q“Qj^ +q 
—i{Fk+k 0 - q — F k+kq _ q — Gk+kq-q + G k + ko - q ){F q — Ft — G g + G g )Q q Qk 0+q 
~i{J~k+ko-q — F k+ko _ q — Gk+ko-q + 5fc+feo-g)(-^ r fco+g — F kq _y q + Gko+q — 5fco+g)Qk 0 +qQk 0 +q 
-2(Fj - G q )(F k+k0 -q + Gk+k 0 -q)Ql Z QZ+^ + 2 ( F lq + Gk-q)(Fk 0 +q ~ G k 0+ ,)Q^Q^+q 
-2(Ft Q+? - S fc0+9 )(F fe _ 9 + G fc _ g )QfQ^ +q + 2(Ft +fco _ g + Gk+k 0 -q)(F q - G q )Q v k z Q % +q 
+ ( — -Fk+ko-q + -^fc+fco-g + Gk+ko-q ~ Gk+k 0 -q)(F q + F^ — G q — Gg) < 9k_ q Qk 0+q 
"^fco+g ~ Gko+q — Gkq+q ) ( F k — q + F k _ q + Gk-q — Gk-q)Qy +ko _ c{ Q kQ ^ 


F(-Fk 0 +q 

F(Fk 0 +q 


Gko+q Gk 0 +q) ( F k - q + F k _ q + Gk-q G k -q)Qq Qk 0 +tJ 
+ ( — F k +ko-q + + Gk+ko-q ~ Gk+k 0 -q){F q + F^ — Gq — G g )Qj^ o+q ^ . 


o+q 


Appendix C: Analysis of the decay and confluence processes 

Omitting the dipolar interaction, the spectrum of Heisenberg AF has the form ek = Dk — DL> 2 k 3 at k 1 (cf. 
Eq. 09). One leads to the following expressions for the decay and confluence processes, respectively, using this 
spectrum: 


Ck Cq Ck-q 


Ck — £q + £q-k 


D(k — <7 — |k — q|) — 3 DL 2 kq(k — q) 
D(k — g + |q — k|) + 3 DL 2 kq(q - k) 


Dka 2 

~MkFr)- 3DLMk ~' l)<0 - 

7. . , + 3DL 2 kq{q - k) > 0, 

2 q(q - k) 


(Cl) 

(C2) 


where q|| and qj. are components of q which are parallel and perpendicular to k (see Fig. [3j , respectively, and we 
assume that q and k are nearly parallel each other (i.e., q» « q and q± Cg). It is seen from Eqs. (IC1I) and (IC2I) that 
both confluence and decay processes are impossible without dipolar forces. 
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1. Decay processes 


Among 8 allowed decay processes (l32l) only the following ones appear to be possible if we take into account the 
dipolar forces: 


4 - e q - e k-q = -3DL 2 k [q(l-q) + C 


£ k - £ q - £ k-q 


f+ — c+ — c~ 

fc k fc q fc k-q 


= -3 DL 2 k 3 g(l - q) + C 


where q = q/k, 0 < q < 1 , and 


= -3 DL 2 k 3 (q{l - q) + C 


A 2 


■ l -g(l-g) 

9(1 - 9 ) 

■ 1 - 9 ( 1 - 9 ) 

9(1 - 9 ) 

■ 1 - 9(1 ~ 9 ) 

9(1 - 9 ) 




m 


2kq(l - q) 


Dqi 


2kq(l - q) 


= 0, 


= 0, 


- r?(l - g) - 


Dq\ 


2kq{l - q) 


= 0, 


c = 


v = 


A 2 

6 L 2 D 2 k 4 ~ 72L 2 J 2 S 2 k i ’ 

woii(v ? k)sin 2 6> k 


72 L 2 Jk 2 

It is seen that Eq. (IC3I) can have a solution if the following inequality holds: 

- (z 2 + C(1 -z)- T)z) < 0, 
z 

where z = q( 1 — q) and 0 < z < 1/4. Solving the quadratic equation, one finds that Eq. (1C8I) is satisfied when 

< q(l -q)< z + , 


At = \ (c + V ± \/(C + vf - 4C ) • 

It is convenient to discuss a limiting case of 77 C that reads 

T A 


k > 


oj 0 D sin 0 k 


(C3) 

(C4) 

(C5) 

(C 6 ) 

(C7) 

(C 8 ) 

(C9) 

(CIO) 

(Cll) 


The opposite limit of 77 <C C has no meaning because it could be realized for £ > 4 only in which case z > 1. One has 
from Eq. (IC10D at ?? C 


2+ « 17, 


C 


(C12) 


The requirement z+ Z- reads 


sin 0 k » 


A 

SW 


(C13) 


It is seen from Eq. (IC7I) that 77 <C 1 if IjCllI) holds. As a result there are two intervals for q inside which inequality 
is satisfied: 


9 e 
9 e 



(C14) 

(CIS) 


It is easy to demonstrate that solutions of Eqs. (IC4I) and (IC5I) exist only for q lying inside intervals (IC'151) and 
(IC14I) . respectively. 
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2. Confluence processes 


Possible confluence processes have the form 

e k+ e q-k~ e q = 3 DL 2 k 3 (q(q - 1) + 

e k + e q -k - < = 3 DL 2 k 3 (q{q - 1) + ^ 1 

e k + e q ~ e q-k = 3DL 2 k 3 ^ q(q - 1) + 


V + 


Dq\ 


2kq(q - 1) 


= 0, 


™' + mh)= 0 ' 


(C16) 

(C17) 

(CIS) 


and there are also other three processes which differ from the presented ones by replacement of q by k — q. Comparing 
Eqs. (IC3I) (IC5I) and (IC16I) (1C18I) one concludes that it is necessary to analyze similar inequality onz = q{q — 1) 


Inequality (IC19I) is satisfied if 


- (z 2 + C(1 + z)~ rjz) < 0. 
z 

z- < q(q~ 1) < z+, 

z± = ^(v-C± 'Jiv- C ) 2 - ■ 


One has in the limiting case of 77 C 

q€ (l + ^l + vj 


(C19) 


(C20) 

(C21) 


(C22) 


so that q ~ k inside this interval. 

The opposite limiting case of 77 -C C is also possible for Eqs. (1C17I) and (jC18l) . This limiting case corresponds to 
q^> k and 


a nr a 

n < \ — n . 7, ■ 

U \ loq L> sin t/k 


(C23) 


In order Eqs. (1C17I) and (1C18I1 have solutions, q should lie between roots of the equation z 2 — qz + ( = 0, i.e., in the 
interval 


qe 



(C24) 


Quantities A and B defined by Eqs. (03 and (PI) . respectively, are related to ( and 77 as follows: A = 77/2 and 
B = 4£/?7 2 . 


Appendix D: Competing single-ion anisotropy and effect of impurities 

We discuss in this appendix the effect of a cubic magnetocrystalline anisotropy on the properties of Heisenberg AF 
with dipolar forces considered in the main text. Assuming for simplicity that S is large, one can model the effect of 
the cubic anisotropy by the following single-ion interaction: 

^ E ((s?) 2 (s ?) 2 + (sf ) 2 {sd 2 + (smsn 2 ). (di) 

i 

We imply below that Uq/SJ ~ |A| <C wo <C J. Let us assume also that the anisotropy constant A' and the spin value 
S' differ from A and S, respectively, at some randomly distributed sites, which concentration is equal to x <C 1. As a 
result, bilinear part of the Hamiltonian (0 acquires the following correction if the staggered magnetization is directed 
along a cube edge: 

W? is = 2AS E a t a k + 2 M' - A)Sj^ X (n)4a n , (D2) 

k 


n 
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where a = S'/S, x( n ) = 1 at sites occupied by impurities, and x(n) = 0 at other sites. Discussion of the Hamiltonian 
"H 2 + 'H. 2 nls can be carried out in the first order in x using the T-matrix approach as it is done, e.g., in Ref. 30 '. The 
situation here is simplified greatly by two circumstances: i) |A|, \A!\ <C J, and ii) sums over k of Green’s functions 
F(w, k), Ft(w,k), t/(w,k), C/(w,k), ,F(w,k), and J^(w,k) at w <C SJ are of the order of ojq/SJ 2 whereas such sum 
for G(w, k) = G(—w, — k)* is much greater being of the order of 1/5 J. As a consequence, the greatest contributions 
from Eq. m arises in £k and £k- Then, one obtains from Eqs. m and m for the correction to 0(w, k) at k <C 1: 




-F 


( w = e ok + *8, k) 


2-\/dk 


= 45 2 Jo A + x{aA ' - A) + 2 x(aA' - A) 2 S— 


SJ 0 


N 


( e 0k + *^) 2 — S 2 (Jq — J 2 ) 


(D3) 


To derive the last term in Eq. m, we set wo = 0 in )>/) q (G(w, q) + G(w, q)) and use Eq. (IA2D for the normal Green’s 
function. The imaginary part of the last term in Eq. (0 determines the magnon damping due to the scattering 
on impurities whereas its real part is negligibly small compared to the second term because |A|,|A'| -C J. Using 
Eqs. (E5l) and (ID3I) . we obtain for the square of the gap in the spectrum 

A 2 = A 2 + 45 2 J 0 (A + x(aA' - A )), (D4) 


where A is the contribution to the gap from dipolar forces given by Eq. ®. The magnon damping due to the 
scattering on impurities is estimated from Eqs. (PI and imi as 


r k 


xS 


(aA' - A) 2 
J 


(D5) 


It is seen from Eq. (ID4I) that the gap in the spectrum can vanish if the anisotropies of dipolar origin (E3l) (accom¬ 
panied by the gap A), A and A! compete. The change of the sign of A 2 signifies that the easy direction switches from 
a cube edge to a cube space diagonal. For instance, the competition arises in TlMni_a;Co x F 3 and RbMni_ x Coa,F 3 
between A < 0 (favoring the cube space diagonals) on the one hand and A! > 0 and the dipolar anisotropy (favoring 
the cube edges) on the other hand. The gap (ID4I) vanishes in this case when x is equal to 


A 2 +45 2 J 0 A 
45 2 Jq(A — a A')' 


(D6) 


Notice that a;* given by Eq. (lD6l is much smaller than unity in TlMni-^Co^Fs and RbMni_ x Coa,F 3 because A' | A\ 
and A 2 ~ 45 2 Jo|A|. 

The above consideration is justified when ej. Tk- It is seen from Eq. (ID5I) that this condition can be invalid for 
some k. For instance, £k S> Tk at x ~ a;* if 


k 


» x 


{aA' - A) 2 
J 2 


(D7) 


The invalidity of the inequality £k ^ Tk can signify a localization of states near the magnon band bottom (see, e.g., 
Ref^£ and references therein). Then, Eq. (ID6I) is just an estimation of the concentration value at which the easy 
direction switches because Eq. (ID4I) is invalid when A is smaller than Tk given by Eq. (ID5I) . In TlMni-^Co^Fs 
and RbMni.^-Co^Fs at x ~ a:*, inequality (ID7I) reads as k 10~ 6 A” 1 . As soon as the great damping due to 
magnon interaction is expected in these compounds at k < 10” 3 4- 10 -4 A -1 , these materials can be suitable for the 
experimental observation of the magnon breakdown discussed in the main text. 

It should be noted also that as soon as defects change the bare spectrum considerably at x « x t and k <C 1, the 
results obtained above must be used with caution. In particular, one cannot exclude the possibility of great spectrum 
change by terms of higher orders in x. It is difficult to analyze the whole series in x but we notice that all the expected 
contributions are small as compared to those considered above due to the smallness of ratios A/J and A'/J. 
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